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We obtain the Na I A8183,8195 absorption line radial velocity curves for the polars ST 
LMi and MR Ser, from which we find the semi-amplitudes to be K a b s —329±6 kms" 1 
and K a b s =289±9 kms -1 respectively. We find that for both systems the effects on 
the Na I absorption lines due to X-rays heating the inner face of the secondary are 
C*~) ' negligible, and so the values obtained for K a i, s can be taken as the true semi-amplitude 

of the secondary star. We then determine the projected rotational velocities, V ro t sini, 
to be 104 ±9 kms -1 and 66±13 kms" 1 for ST LMi and MR Ser respectively which 
enables their mass ratios to be calculated. For ST LMi and MR Ser we find the mass 
ratio to be 0.22±0.04 and 0.10±0.05 respectively; values which are significantly dif- 
ferent only at the 94 percent level. Using the value for the orbital inclination derived 
from polarimetric measurements, we determine the mass, and the orbital and rota- 
\Q ■ tional velocities of the secondary stars. These are significantly different at less than 

f^*) ' the 90 percent level. However, if the limb darkening is the same in both objects, these 

^sO . quantities are significantly different at the 96 percent level. 

<3\ ' We present Doppler maps of the Na I absorption and Ca II emission in ST LMi 

"pH \ and MR Ser. In both systems the Na I absorption covers the secondary star. In ST 

O LMi the Doppler image of the Ca II emission shows that it originates at the inner face 

I 1 of the secondary star. In MR Ser, however, the emission lies close to but not on the 

secondary star. 

CO 



We show that "spike" in the orbital period distribution of polars is a significant 
feature, although the discovery of only one more system with a period outside the 
"spike" would decrease its significance below a 99 percent confidence level. We conclude 
that, even if the limb darkening coefficients for the secondary stars in ST LMi and MR 
Ser are the same, we cannot rule out the two systems having identical parameters. 
Therefore our observations are compatible with the theory explaining the "spike" in 



C3 ' the period distribution of the AM Hers. 
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1 INTRODUCTION 

AM Her type systems consist of a magnetic white dwarf 
primary star accreting material from a late type M dwarf 
secondary star. The white dwarf primary star has a strong 
enough magnetic field (B>10 MG) to lock it into syn- 
chronous rotation with the orbit. If one compares the or- 
bital period distribution of the AM Hers with that of the 
non-magnetic CVs, then one finds that the main features 
of the CV distribution re-appear in the form of a minimum 
period P or (,=78 m, and the period "gap" between P or b~2-3 
hrs (Wheatley 1995). However, there are three obvious fea- 
tures in the AM Her period distribution (1) a concentration 
of systems near P or i,=114 m (the "spike"), no such feature 



occurs in the distribution of weakly magnetic CVs; (2) a 
long period cut-off, significantly shorter than that for CVs, 
and (3) more systems in the period "gap" than for CVs. 

Cataclysmic variables evolve from longer to shorter pe- 
riods, ceasing mass transfer as they pass through the period 
"gap". The period "gap" has been modelled in terms of an 
interruption of orbital angular momentum loss at 3 hrs, al- 
lowing the secondary star to relax into thermal equilibrium 
and detach from its Roche lobe. Systems regain contact at 
2 hrs after orbital angular momentum loss through gravita- 
tional radiation. The "spike" in the polar distribution has 
been explained as being the period at which most polars 
come out of the period "gap" (Hameury, King & Lasota 
1990), However, it has been shown that for the "spike" to 
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be so narrow, for those systems in the "spike" , the masses of 
the secondary stars must ALL be 0.200 ± 0.002 M© AND 
the masses of the white dwarfs must only differ by ~ 0.05 
M Q (Ritter & Kolb 1992). Therefore not only do we ex- 
pect the mass ratios for the systems in the "spike" to be 
approximately the same, but the general system parameters 
such as the orbital velocities should also be almost identi- 
cal. To test this theory the system parameters of polars and 
in particular those in the spike must be determined more 
accurately than before. This can be done by measuring the 
radial velocity of the secondary star, K2, and its rotational 
velocity, V ro tsini. In the red region of the spectrum, where 
the M dwarf is expected to contribute significantly to the 
observed flux, the spectra of polars are usually dominated 
by cyclotron radiation, which originates in the post-shock 
flow near the magnetic pole of the white dwarf, thus mak- 
ing observations of the red component difficult. However, 
when they enter the low state and the accretion deceases, or 
in parts of the orbit where the accreting region is occulted 
it becomes much easier to observe the secondary star (see 
Cropper 1990 for a review). 

In this paper we report on the use of medium resolution, 
phase resolved spectroscopy to study the orbital motion of 
the secondary stars in ST LMi, MR Ser, AM Her and 1H 
1752+081. We measure V rot sini and K 2 for MR Ser and ST 
LMi and deduce their system parameters. We also measure 
Vrot sin i for AM Her and attempt to do so for 1H 1752+081. 
The polars MR Ser, ST LMi and 1H 1752+081 lie in the 
"spike" . We also discuss the likelihood of the "spike" being 
a genuine feature and if our observations of ST LMi and 
MR Ser are compatible with theoretical explanations of the 
"spike". 



2 OBSERVATIONS AND DATA REDUCTION 

MR Ser, ST LMi, AM Her and 1H 1752+081 were observed 
with the ISIS spectrograph attached to the 4.2-m William 
Herschel Telescope on the nights of 1994 March 24/25 and 
25/26. Here we only present the results from the red arm. 
The EEV CCD was used with the R600R grating giving a 
dispersion of 0.74 A pixel -1 . On the first night the wave- 
length coverage was 7800-8586A but on the second night 
the central position was moved to give 8030-8888A cover- 
age. Our principal aim was to determine radial velocities 
and to measure the rotational broadening of the absorption 
line, so we used a narrow slit (1.0 arcsec) which resulted in 
a spectral resolution of 1.9A (FWHM) =70 kms -1 . Further 
observational details are contained in Table 1. The seeing 
on the first night was steady at about 1 arcsec, and on the 
second night it was typically between 0.5 and 1.2 arcsecs. 
The spectral type of the secondary star in AM Her has been 
determined by many authors to be M5 (see Southwell et al. 
1995 and references within). Mukai & Charles (1987) deter- 
mined the spectral type of the secondary stars in ST LMi 
and MR Ser to be M5.5. We therefore also observed Gl 406 
(M6V) and Gl 473 (M5V) which we use for cross-correlation 
in the next section. We took care to observe this star through 
the slit of the same width as used for the two targets. The 
flux standard Feige 34 (Oke & Gunn 1983) was observed 
using wide and narrow slits on each night to correct for the 
instrumental response. Cu-Ar arc exposures were taken ev- 



ery 30-40 min in order to calibrate the wavelength scale. We 
binned the CCD pixels on readout into groups of three pixels 
along the slit. The effect of this was to reduce the readout 
noise and deadtime of the CCD. 

Sensitivity variations were removed with a balance 
frame prepared from tungsten-lamp flat-fields. Third order 
polynomial fits to the sky were subtracted, and raw spectra 
were then extracted using the optimal algorithm of Hornc 
(1986), which weights the pixels along the spatial direction 
in order to obtain the maximum signal-to-noise ratio. The 
arc spectra were extracted with the same weights and at 
the same position as the target, and fitted with third-order 
polynomials. The rms scatter of the fits was ~ l/20th of a 
pixel. With four coefficients fitted to about 20 lines in each 
case, the statistical uncertainly in the wavelength scale is 
therefore of the order l/45th of a pixel. This corresponds to 
an uncertainty of 0.6 kms -1 . 

Atmospheric features were removed from the red arm 
spectra by using a wavelength calibrated water template star 
produced from the observed standard star. The template 
star was then corrected to the zenith distance (z) of the 
object by multiplying by a factor of (sec 2) ' 6 . A comparison 
of the tabulated absolute flux values with a spline fit to the 
continuum of the flux standard provided a correction for the 
large-scale instrumental response. 



3 RESULTS 

3.1 Radial velocities of ST LMi and MR Ser 

We obtained absorption line radial velocity curves using the 
Na I A8183,8194 doublet and also emission line radial ve- 
locity curves using the Ca II A8498, 8542 lines. K abs and 
K erni are defined as the semi-amplitude of the absorption 
and emission line radial velocity curves respectively. The he- 
liocentric radial velocities of MR Ser and ST LMi are shown 
in Figures 1 and 2. 

In order to measure the absorption line radial velocity 
curves we cross-correlated the individual spectra with the 
template star spectrum (Tonry & Davis 1979), only using 
the wavelength range close to the Na I absorption lines. 
Prior to cross-correlation, the spectra were interpolated onto 
a logarithmic wavelength scale. Three spectra from ST LMi 
and eighteen from MR Ser were too noisy to obtain a reliable 
cross-correlation peak, and were not used for the fit. Gl 406 
(M6V) was used as the template spectrum. We find K abs to 
be 329 ± 6 km s -1 and 289 ± 9 km s" 1 for ST LMi and MR 
Ser, respectively. 

For each of MR Ser and ST LMi, to obtain the tem- 
plate spectrum for the Ca II triplet emission lines we ob- 
tained a preliminary template by averaging the spectra using 
K abs . The summed spectrum showed some double peaked 
structure implying that the velocity used to correct for the 
Doppler motion was incorrect, and was not the velocity of 
the regions producing the Ca II emission. We then per- 
formed the cross correlation to obtain K ern i, only using the 
wavelength range which contained the Ca II A8498, 8542 
emission lines. This value was then used to Doppler cor- 
rect the original spectra to give the final template spectrum 
for the Ca II emission, in which the emission was narrow 
and single peaked. We find K emi to be 213 ± 4 km s _1 and 
173 ± 1 km s _1 for ST LMi and MR Ser, respectively. 
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Because cataclysmic variables have very short periods, 
it is generally believed that any initial non-circularity would 
have been quickly removed by tidal forces between the red 
star and the white dwarf, and that the present orbits are in- 
deed circular. However, Davey and Smith (1992) argue that 
the radial velocity curves may still be distorted from a pure 
sine wave by geometrical distortion and heating of the sec- 
ondary star by its companion, causing the centre of light 
given by the strength of the Na I doublet to differ from the 
centre of mass. The effects of this can be represented by 
allowing for a phase shift in the sine curve, or more gener- 
ally by introducing a fictitious eccentricity. They describe a 
procedure for detecting any effects of heating on the radial 
velocity curve. Firstly one must check for the significance of 
a fit with an eccentric orbit. If the fit is not significantly bet- 
ter than a purely sinusoidal fit, then the semi-amplitude of 
the curve K a b s measured from the absorption features rep- 
resents a measure of the true semi-amplitude of the radial 
velocity curve of the secondary star Ki- 

For both systems we found that the fit to the absorp- 
tion line radial velocity curve with an eccentric orbit was not 
significantly better than that with a circular orbit. We ob- 
tained eccentricities of 0.027±0.017 and 0.075±0.046 for ST 
LMi and MR Ser, for which the significance of the eccen- 
tricity parameter was less than 50 percent. Therefore the 
measured value of K a b s from the absorption lines can be 
taken to represent the true value of Ki. After adding the 
heliocentric radial velocity of the template star (13 kms -1 ; 
Gliese 1969) we determined Ki using a circular fit of the 
form 

V(<t>) = 7 + K 2 sin(27T0) (1) 

to the absorption and emission line heliocentric radial ve- 
locities (see Table 2). Here 7 is the systemic velocity and <f> 
is the phase angle of the observations with 0=0.0 at infe- 
rior conjunction (i.e. when the red dwarf is in front of the 
white dwarf) . The ephemerides derived for ST LMi and MR 
Ser (see Table 2) agree well with those of Mukai & Charles 
(1987). Figures 3 and 4 show the Doppler corrected average 
spectra of the two targets. 

As another check to see if there is any effect due to 
heating, we fitted the Na I radial velocities measured only 
between phases —0.2 and +0.2, when the secondary star is 
viewed from the back and its appearance should be rela- 
tively unaffected by the effects of irradiation. We find that 
the value for K abs does not change within the uncertainties. 
Our derived values for K2 from the sinusoidal fit to the radial 
velocity curve are comparable with those derived by Mukai 
& Charles (1987). However there seems to be some discrep- 
ancy in the value for 7, which may be explained by the fact 
that their are some inconsistencies between the published 
radial velocity for most red dwarfs (Friend et al. 1990). In 
order to check our 7 velocity estimates we performed the 
cross correlation analysis using a fake template star, hav- 
ing the instrumental resolution of the data, and the main 
absorption and emission lines present at their rest wave- 
lengths. This implies that the fake template star has no 7 
velocity associated with it. We found that our derived pa- 
rameters agreed well with those obtained earlier. Despite 
this, it should be noted for the emission line case, that using 
the Doppler averaged sum of the target star as a template 
star for the emission line data may in some way bias the de- 



termination of the emission line 7 velocity. The result may 
depend on the phase chosen as the reference phase for the 
average spectrum. Therefore, the only reliable estimate of 
the 7 velocity from the emission lines is that obtained using 
the fake template star. 

3.2 AM Her and 1H 1752+081 

Our data were less extensive for AM Her and 1H 1752+081. 
For AM Her we cross-correlated the template star Gl 473 
(M5V) with the target spectra in order to determine the 
velocity shift. We then Doppler corrected and summed the 
spectra. The mean spectrum is shown in Figure 5. For 1H 
1752+081 the signal-to-noise ratio of the data was such 
that we were unable to Doppler correct the spectra to the 
Na I radial velocity. We therefore merely show the summed 
spectrum in Figure 6. Only low resolution blue spectra exist 
for this source (Ferrario et al. 1995). Many emission lines can 
be seen in both objects such as He II A8246, Ca II A8498, 
8542, 8662 and the Paschen lines A8598, 8750. In AM Her 
the Na I absorption lines are strong whereas they are weak 
in 1H 1752+081. Na I has not been observed in 1H 1752+081 
previously. 

3.3 The rotational broadening and mass ratio 

In order to determine the rotational broadening, V ro t sin i, 
we followed the procedure described by Marsh, Robinson 
and Wood (1994). We Doppler corrected and summed the 
spectra for our individual targets using K abs . We then fit- 
ted the excess light from the system i.e. that not arising 
from the secondary star, plus a constant times a rotation- 
ally broadened version of the template star, to the object's 
spectrum. The template star was broadened by convolution 
with the rotational profile of Gray (1976) assuming a lin- 
earized limb darkening coefficient of u = 0.72, taken from 
Al-Naimiy (1978). The broadened templates and the excess 
light are shown for the individual objects in Figures 3, 4 
and 5. A high pass filter was then applied to the difference 
spectrum to remove any large-scale spectral differences and 
then the \ 2 was computed. Fitting a parabola to the distri- 
bution of x 2 versus broadening values, we obtain a minimum 
X 2 for each target and so the best value for the rotational 
broadening. 

For ST LMi and MR Ser we obtained V ro t sin i=105±7 
and 66+12 km s _1 respectively, and for AM Her we obtained 
l^ot sin i=70±ll kms -1 . The latter is consistent with that 
obtained by Southwell et al. (1995). We were not able to de- 
termine Vrotsmi in 1H 1752+081. The 1-a uncertainties in 
Vrot sin i were derived by forcing the minimum x 2 to increase 
by 1 (Lampton, Margon & Bowyer 1976). 

The analysis above gives the uncertainties assuming the 
limb darkening is known. However, in reality this is not the 
case. We therefore repeated the above procedure using no 
and full limb darkening. For u = 0.0 we find V ro t sin i — 
98 ± 7, 60 ± 12 and 65 ± 10 kms _1 for ST LMi, MR Ser 
and AM Her respectively. For u = 1.0 we find V ro t sin i — 
110 + 7, 71 + 12 and 72+12 kms -1 for ST LMi, MR Ser and 
AM Her respectively. We therefore adopt V ro t sin i— 104+9 
km s" 1 and V TOt sini=66±13 km s" 1 for ST LMi and MR Ser 
(see Figures 3 and 4). For AM Her we find V rot sini=68±12 
km s -1 (see Figure 5). 
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Armed with K 2 and V ro t sin i we can now determine 
the mass ratio q(— M2/M1) Assuming that the companion 
is tidally locked, spherically symmetric and fills its Roche 
lobe, the rotational broadening provides a measurement of 
the binary mass ratio through the result (see e.g. Wade & 
Home 1988) 



Vrot sini = 0AQ2K 2 q ir '\l + q) 



n2/3 



(2) 



Substituting our values for Vrot sin i and K2 in equation 
(2) we find the mass ratios of ST LMi and MR Ser to be 
0.22±0.04 and 0.10±0.05 respectively. (Since our results for 
AM Her are consistent with those of Southwell et al. (1995) 
and they have already determined the system parameters we 
do not do so here). The correction of V ro t sini and q due to 
the Roche geometry only introduces an error of a few per- 
cent (Marsh et al 1994; Southwell et al 1995). We performed 
a t-test of the values obtained for the mass ratio for ST LMi 
and MR Ser and find that the values for the mass ratios are 
only significantly different at the 94 percent level (1.9<r). 

However, if we assume that the two systems have the 
same limb darkening coefficient, for example u — 0.72, we 
then find the uncertainties in q for ST LMi and MR Ser are 
reduced to 0.03 and 0.05 respectively. Using these uncer- 
tainties, we then find that the mass ratios are significantly 
different at the 96 percent level. We thus conclude that even 
if we assume that the systems have similar limb darkening 
we cannot show at the 99 percent confidence level that the 
two systems have different mass ratios. 



3.4 Doppler images 

The emission lines from cataclysmic variables are broadened 
by Doppler shifting from the accretion disc (in weakly mag- 
netic systems), component stars and mass transfer. Marsh 
& Home (1988) have shown that it is possible to derive the 
distribution of emission from observations of the profile cov- 
ering the binary orbit. Here we use the same technique but 
the features we are interested in are narrow and arise from 
the companion star or the accretion stream. The images are 
presented in velocity space rather than position coordinates, 
since there is no unique transform between velocity and po- 
sition in accreting binaries. 

Figures 7 and 8 show the Doppler maps of the 
Na I absorption and Ca II emission distribution for ST LMi 
and MR Ser respectively. We also place on the maps the 
path of the gas stream and the shape of the secondary star's 
Roche-lobe for the given values of K2 and q. The Na I maps 
for both systems show absorption covering most of the sec- 
ondary star, and the Ca II emission arising from regions 
near the inner Lagrangian point (£1). It should be noted 
that the resolution of the Doppler maps is limited by the 
velocity dispersion of the data (see section 2) . 

The emission-line region is displaced towards the centre 
of mass of the binary system, making the observed veloc- 
ity amplitude smaller than the true secondary star velocity. 
This "K2" correction (Wade & Home 1988) can be expressed 



K2 = 



(3) 



the secondary star, and < f < 1 expresses the distribu- 
tion of line emission over the surface of the secondary star. 
Placing the emission entirely at the L\ point i.e. f=l, for 
ST LMi using K emi =213± 4 kms" 1 we find K 2 =313±13 
kms -1 . This agrees with the value for K2 obtained from 
the absorption lines. For MR Ser, the "K2" correction gives 
a value for K2 much smaller that that obtained from the 
absorption line radial velocity curve. One requires q >0.35 
in order to obtain the observed K2- These results are con- 
sistent with the fact that the emission in the Doppler maps 
in ST LMi lies on the secondary star near L\, while in MR 
Ser it is close to L\ but not on the secondary star. 

The accretion stream can be divided into two distinct 
parts. The first part is free-falling in the orbital plane and 
the second is where the stream grabs onto the magnetic field 
lines at the magnetospheric radius. The phasing of the first 
narrow emission component would be offset from that of the 
absorption lines arising from the secondary star. The mag- 
netic field in all AM Her type systems in which it has been 
measured has been found to be similar, including ST LMi 
(Schmidt, Stockman & Grandi 1983) and MR Ser (Wickra- 
masinghe et al. 1991). For systems with periods less than 
2 hrs like ST LMi and MR Ser the magnetospheric radius 
is very likely to be close to the L\ point. Observations of 
ST LMi show that the white dwarf magnetospheric radius 
is 0.79-0.93i?Li and in MR Ser R m a g is expected to be very 
similar (Mukai 1988; where Rli is the distance of the in- 
ner Lagrangian point from the white dwarf). In MR Ser the 
emission we observe in velocity space is very close to the 
Li point. It could arise from an extended emission region 
around the L\ point. If it is due to emission from the region 
where the free-falling gas stream joins onto the magnetic 
field lines of the white dwarf, then one would expect it to lie 
on the gas stream in the Doppler map. This is not the case. 



4 DISCUSSION 

4.1 The mass of the binary components 

Given our values of K2 and q combined with the orbital pe- 
riod, P, and the binary inclination, i, we can also determine 
the masses of the compact object and its companion using 



2ttG 



Mi sin 3 



(4) 



l-f(l + q)(R 2 /a) 
where R2 is the distance from the L\ point to the centre of 



The binary inclination for MR Ser and ST LMi can be ob- 
tained from polarimetric observations. For ST LMi and MR 
Ser we use inclination angles of 56° ±4° (Cropper 1988) and 
45° ± 5° (Brainerd & Lamb 1985) respectively. We find that 
the masses for the secondary stars in ST LMi and MR Ser 
are M 2 =0. 17+0.07 M Q and M 2 =0. 07+0.08 M Q respectively 
(see Table 2) . Using the values for the inclination we can de- 
rive the orbital velocity, V or b = K2 / sin i, and the rotational 
velocity, V ro t, of the secondary star. For ST LMi and MR Ser 
we find Vorb to be 397±20 kms _1 and 409±38 kms -1 , and 
Vrot to be 125+12 kms -1 and 93+20 km s -1 respectively. If 
we assume that the two systems have similar limb darken- 
ing, then the uncertainty in the V ro t measurements reduces 
to 10 km s -1 and 19 km s -1 for ST LMi and MR Ser respec- 
tively. Similarly, for V or t the uncertainties become 20 and 
38 kms _1 and for M 2 0.06 and 0.08 M Q . The parameters 
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are still not significantly different at the 90 percent confi- 
dence level. 

4.2 The "spike" 

The sample of known AM Hers has more than doubled over 
the last 5 years, largely as a result of the ROSAT extreme 
ultraviolet and soft X-ray surveys. To date there is a total 
of 43 systems (see Table 3). In Table 4 we list some details 
of how the AM Hers are distributed. We set the lower and 
upper edges of the period "gap" to lie at 126 and 192 mins 
respectively. The upper edge of the whole distribution is 
taken to be at 300 mins. The outlier RXJ05I5+0I with an 
orbital period of 479 mins is excluded from the analysis that 
follows because of its unusual nature. Its low magnetic field 
strength is somewhat of a surprise, because for a given white 
dwarf mass one would expect that the field strength required 
to enforce synchronism should increase significantly as the 
orbital period increases (Buckley & Shatter 1995). 

There are two questions we want to answer in order to 
explain the observed characteristics of the period distribu- 
tion (see Figure 9). 

• What is the probability of the observed number of sys- 
tems lying in the period "gap"? 

• What is the probability of the "spike" occurring by a 
statistical fluke? 

There are 42 AM Her systems in the period range used 
(see above) , of which 4 are in the period "gap" which lies 
between 126 and 192 mins. Following the calculation of 
King 1994 (private communication) if the AM Hers are dis- 
tributed uniformly in the range 78-300 mins, the probability 
of finding at least 4 systems in the period "gap" is given by 

P(4 in "gap") = N C 4P N - 4 (1 - p) 4 (5) 

where N C± is the number of ways of choosing 4 objects from 
a set of N, N is the total number of systems in the period 
range and p is the probability of an AM Her avoiding the 
"gap". Using iV=42, we find p=(l-66/222)=0.703, which 
implies P(4 in "gap") = 1.3xl0 , i.e. the observed number 
of systems in the gap is inconsistent with the distribution 
being uniform at the 99.9 percent level. 

The pre-ROSAT data contained 12 members below the 
period "gap", of which 6 were in a 2 min interval (the 
"spike"). In the current sample the number of "spike" mem- 
bers has risen to 9 at the cost of the "spike" broadening to 3 
min. Suppose that there are M AM Her systems below the 
period gap, uniformly distributed over the period range of 
R mins, so that the probability of any one system being in 
a given interval of 3 mins is 3/7?. We want to know what is 
the probability of the remainder M-l being distributed in 
such a way that r — 1 of them (a total of r) will lie in the 
same 3 min interval, for any initial choice of interval. One 
has to add these results for each of the initial choices, which 
means multiplying by M. Thus 

P{«spike") = M M ~ 1 C r -i p r - 1 {l-p) M - r . (6) 

With 7?=38 mins corresponding to a range 78-116 mins 
(the most conservative), M=25 and r~9 systems in the 
"spike", we find P( "spike" )=0. 0068, i.e. the probability that 



the spike is genuine is 99.3 percent. Using the range 78- 
126 mins, M—27 and r=9 systems in the spike we find 
P( "spike" )=0. 0030 with a probability of 99.7 percent that 
it is genuine. The old (pre ROSAT launch) sample had 12 
systems with 6 in a 1.2 min interval. Using data before the 
ROSAT detections, we find the probability of the "spike" 
occurring randomly was 0.00014. 

Although the significance of the "spike" has reduced 
since the ROSAT mission, the "spike" is still a significant 
feature of the period distribution of the AM Hers. However, 
it would take only one new system in the period range 78- 
116 mins but outside the "spike" to raise the probability 
that a 9 member "spike" could be a random event to more 
than 1 percent. 

It should be noted that equation (6) is the probability 
that any "spike" with an interval of 3 mins occurs below 
the "gap" . The question we should next ask, is what is the 
probability of the "spike" with r members occurring in a 
particular 3- min interval at 113-116 mins. This is given by 

P( u spike" at 113 - 116 mins) = M C r p r {I -p) M ~ r (7) 

Using R=38 mins (the most conservative), M—25 and 
r=9 systems in the "spike", we find P( "spike" at 113-116 
mins)=0.00007, i.e. 99.99 percent. Thus the "spike" in the 
period distribution occurring at 113-116 mins is a significant 
feature. 

If we compare the secondary star masses derived earlier 
for ST LMi and MR Ser, with that predicted by the "spike" 
theory of 0.200 ± 0.002M Q (Ritter and Kolb 1992), we find 
that we cannot show with certainty, i.e. at the 99 percent 
level, that the two systems have masses different from this 
theoretical mass. If we accept the spike theory, then the two 
systems should have almost identical system parameters, in- 
cluding q, M2, Vorb and V ro t- The values we have obtained for 
all the system parameters for ST LMi and MR Ser are such 
that we cannot rule out that the two systems have the same 
parameters. Even if the limb darkening coefficients are the 
same q can only be shown to be different at the 96 percent 
level. Our results are therefore compatible with the theoreti- 
cal explanation of the spike given by Ritter and Kolb (1992). 
Given a value for the limb darkening the dominant source of 
uncertainty in the deduced system parameters, other than 
q, is the orbital inclination, i. To improve the sensitivity of 
this test of the spike theory a more accurate determination 
of i is required. Very high resolution spectroscopy is needed 
to reduce the error in V ro tsmi. More systems in the spike 
should also be investigated. 



5 CONCLUSIONS 

We have obtained data of medium spectral and high time 
resolution covering the orbital cycles of the polars ST 
LMi and MR Ser. Our spectra allow us to determine the 
dwarf star's radial velocity semi-amplitude, K abs , using the 
Na I A8183,8195 absorption features. We find that for both 
systems the effects of irradiation on the Na I lines are neg- 
ligible, since there is no distortion of the radial velocity 
curves from a pure sine wave. Therefore K a bs can be taken 
as the true semi-amplitude of the secondary star K2. We 
find that for ST LMi K abs = 329 ± 6 km s" 1 and for 
MR Ser K aba = 289 i 9 km s _1 . Doppler maps of the 
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Na I distribution for the two systems are very similar, cov- 
ering most of the secondary star. However, in contrast, the 
Ca II maps are somewhat different. In ST LMi, the map 
shows the emission arising from the secondary star at the 
inner Lagrangian point, probably due to heating of the inner 
face of the secondary star, whereas in MR Ser the emission 
is observed to arise from regions close to but not on the sec- 
ondary star. For ST LMi if we assume all the emission arises 
from the L\ point, then the semi-amplitude of the emission 
line radial velocity curve gives a value for K2 consistent with 
that obtained from the absorption line radial velocity study. 

The spectra of the secondary stars are significantly 
broadened by rotation. We determine this projected rota- 
tional broadening which enables us to determine the mass 
ratios of the systems. For ST LMi and MR Ser we find 
g=0.22±0.04 and 0.10±0.05, respectively. Using the val- 
ues for the inclination derived from polarimetric measure- 
ments, we find M 2 =0.17±0.07 M©, V orh = 397 ± 20 km s _1 
and Vrot = 125 ± 12 km s" 1 for ST LMi. For MR Ser 
we find M 2 =0.07±0.08 M , V or b = 409 ± 38 km s _1 and 
Vrot = 93 ± 20 km s _1 . We conclude that we cannot rule 
out the two systems having the same system parameters, 
therefore our observations are consistent with the theoreti- 
cal explanation of the spike. 

Finally we show that at present the "spike" is a signifi- 
cant feature of the period distribution of the AM Hers. How- 
ever, the discovery of only one more system with a period 
that does not lie in the "spike" would decrease the signifi- 
cance of the spike below a 99 percent confidence level. 
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Table 1: Journal of observations 



ST LMi MR Ser AM Her 1H 1752+081 

Night of observation 2 1,2 2 2 

Number of spectra 88 80 4 4 

Exposure time 170s = 0.025</> 170s = 0.025</> 280s = 0.025</> 170s = 0.025</> 



Table 2: Derived parameters for ST LMi and MR Ser, Uncertainties are l-cr assuming 
unknown limb darkening 



Parameter 



ST LMi 



MR Ser 



Period (days) 

T (HJD - 2440000) 

7 abs (km s" 1 ) 

7e m i (km s" 1 ) 

Kf s (km s" 1 ) 

Kl mi (km s" 1 ) 

V ro t sin i (km s -1 ) 

Fraction of light from secondary 

q 

f(M!/M ) 



0.0791±0.0003 
9438. 5546±0. 0003 
3±4 
-5±3 
329±6 
213±4 
104±9 
0.76±0.02 
0.22±0.05 
0.29±0.02 



0.0789±0.0001 
9437.696±0.001 
-50±7 
-51±2 
289±9 
173±1 
66±13 
0.20±0.01 
0.10±0.07 
0.20±0.02 



Using the value for the inclination 
derived from polarimetric measurements 



Mx (M ) 
M 2 (M ) 
R 2 (R©) 
a (R ) 
Vorb (km s" 
V rot (km s" 



56±4 
0.76±0.30 
0.17±0.07 
0.20±0.04 
0.75±0.15 
397±20 
125±12 



45±5 
0.67±0.70 
0.07±0.08 
0.15±0.08 
0.70±0.36 
409±38 
93±20 



1 



Table 4: AM Her period distribution parameters 



Number of systems below the period "g 


jap" 


27 


Number of systems above the period "g 


jap" 


11 


Number of systems in the "spike" 




9 


Number of systems in the period "gap" 




4 


Total number of systems 




42 


Period range for "spike" (mins) 




113-116 


Period range for "gap" (mins) 




126-192 


Period range for whole distribution (mins) 


78-300 



1 



Normalised flux (arbritary units) 

o o j- k j- l ro 

b bi b cji b 




Table 3: Currently known AM Her systems in order of increasing orbital period 



TV a in f 


P t ( m ml 
1 orb ^iiiin ^ 


1 ( HP&VPf"^, 


RX0132-65 


78.0 




RE1307+535 


79.7 




RX1015+09 


80.0 




EF Eri 


81.0 


58±12 


RX0153-59 


89.0 




DP Leo 


89.8 


76±10 


RE1149+28 


90.0 




RE 1844-74 


90.0 




RX1 957-57 


99.0 




VV Pup 


100.4 


76±6 


V834 Cam 


101.5 


45±9 


RE0453-42 


102.0 




EP Dra 


104.6 


~ 80 


RE0501-03 


107.0 




RX1002-19 


107.0 




CE Cm 


108.6 




1H 1752+081 


113.0 


75-90 


RE1802+18 


113.0 


<80 


BL Hyi 


113.6 


70±10 


MR Ser 


113.6 


45±5 


ST LMi 


113.9 


56±4 


EK UMa 


114.5 


56±19 


WW Hor 


114.6 




AM UMa 


114.8 


65±5 


1 1 1 3+432 


115.9 




RE2 107-05 


125.0 


80±5 


EU Cnc 


125.4 




UZ For 


126.5 


150±12 


RE0531-46 


133.0 




RE1938-46 


140.0 




Drissen V211R 


160.0 




BY Cam 


199.3 




AM Her 


185.6 


30±5 


V1500 Cve: 


201.0 




RE1940-10 


202.2 




RX0929-2404 

llivUJui/ till/1 


203.4 




RX1007-20 


208.0 




RX2316-05 


209.0 




QQ Vul 


222.5 


60±14 


EXO 0329-26 


228.4 




RX1313-32 


252.0 




RX0203+29 


275.5 




RE0515+01 


479.0 





For systems discovered before the ROSAT mission, the data has been taken from Cropper (1990). For 
the post ROSAT systems, data has been taken from Ritter & Kolb (1995), Beuermann & Thomas (1993), 
Kolb k de Kool (1994). 
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